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Abstract 

This document considers the formation control problem for a group 
of non-holonomic mobile robots under time delayed communications. 
The agents are assumed to be working under a directed and fixed com¬ 
munication topology. A recently developed formation control tech¬ 
nique, based on a consensus protocol for linear systems is revisited. 
The stability analysis of this control technique, in the domain of the 
time delays, is performed using the CTCR paradigm under the SDS 
domain approach. This analysis leads to an exact declaration of the 
stability boundaries in the domain of the delays. This linear control 
law is adapted to the nonlinear dynamics of the non-holonomic carts 
by means of feedback linearization. Simulation results show the appli¬ 
cability of this method to the class of agents considered here. 


1 Introduction 

Cooperation among multiple agents is a research topic which has attracted 
a lot of attention in the past two decades. In order to perform a coopera¬ 
tive task, the members of a team have to reach an agreement on a certain 
variable. If they want to rendezvous at some point, they must agree upon 
the location of that point. If they are to travel together, they must select 
a heading and speed. This gives birth to the consensus problem. Since the 
original agreement protocol of Vickseck et al. [I], expanded by the work 
of Jadbabaie et al. [2], many authors have worked to extend the theory 
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and applications of the consensus problem. Olfati-Saber and Murray [3] set 
up a formal framework to treat the agreement problem, highlighting the 
importance of algebraic graph theory in this topic. They also considered 
the case in which time delays are present in the communication channels. 
This pioneering work was limited to agents modeled by single integrator dy¬ 
namics, and it was soon expanded to consider cases including higher-order 
dynamics HHZj and multiple time delays p{lo]. 

Some researchers have also studied the more general formation genera¬ 
tion or flocking problem mm- Yao et al. m followed an approach based 
on artificial potentials, whereas Olfati-Saber [12] used a consensus term for 
velocity agreement and a distance based repulsion term to create the spacing 
among agents. This methods are heavily non-linear, despite the fact that 
they are applied to agents with simple single or double integrator dynamics. 
Following a different approach, Lin et al. |7| transformed a consensus algo¬ 
rithm into a formation control protocol by including the desired inter-agent 
distances in the control protocol. 

Cepeda-Gomez and Olgac m, also expanded a consensus control law 
PfTO] to transform it into a formation logic. The last two protocols preserve 
the linearity of the control logic, and therefore allow the use of well known 
techniques for the stability analysis with respect to the time delays. All 
these works, however, have considered agents with very simple dynamics. 

Since real life mobile agents have dynamics more complicated than the 
single or double integrator models commonly used, some authors have ex¬ 
panded the field to include the problem of cooperation of non holonomic 
agents. The works of Zhai et al. Listman |15] and Ghommam m 
present different approaches to solve this problem. However, very few works, 
besides that of Dong and Farrell have treated the consensus or forma¬ 
tion control problems for non-holonomic agents under time delayed com¬ 
munications. To study the stability with respect to the delay, they use an 
approach based on the D-subdivision method, and obtain a delay margin. 
This result is valid only if all the eigenvalues of the Laplacian matrix gen¬ 
erated by the communication topology are real. 

In this work, we depart from a consensus protocol affected by two time 
delays, expand it to include a constant forcing term to create the formation 
and then deploy it into non-holonomic agents, more specifically, Differential- 
Drive Mobile Robots, DDMR. A feedback linearization procedure allows us 
to deploy the original double integrator control protocol to this class of 
agents, and the Gluster treatment of Gharacteristic Roots Paradigm [TSlfT^ 
allows us to find non-conservative stability bounds in the domain of the time 
delays. 
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To summarize the contribution of the present paper is the following. Us¬ 
ing a feedback linearization model, we adapt a recently developed formation 
control technique to a set of non-holonomic robots. The formation con¬ 
trol technique is based on a consensus protocol and assumes two different 
time delays in the sensing of positions and velocities. By using the Cluster 
Treatment of Characteristic Roots methodology, we determine, exactly and 
exhaustively, the delay combinations that guarantee stability, i.e., conver¬ 
gence to formation in a finite time. To the best of our knowledge, this is the 
first time such formation problem under time delays is solved in an exactly 
manner for non-holonomic agents. 

The rest of the paper includes the description of the consensus protocol 
in Section [21 the stability analysis with respect to the delays, using CTCR, 
in Section [3l the formation control algorithm in Section HI the feedback 
linearization procedure in Section [5] and some simulation examples in Sec¬ 
tion m Conclusions and future work directions are presented in Section 0 
Throughout the text boldface notation is used for matrices (uppercase) and 
vectors (lowercase), whereas italic symbols represent scalar quantities. 

2 Consensus Protocol 

We start with the consensus protocol introduced in m- It considers a group 
of n agents with double integrator dynamics of the form Xj (f) = Uj (t ), where 
Xj (t) E represents the position vector of the f-th agent and Uj (t) E M™' 
is the input vector. It is assumed that each one of the position coordinates 
is independent of the other components of the x vector, and that the Ui^k 
input affects only the dynamics of the Xi^k position, i.e., the dynamics of all 
the position coordinates are decoupled. Taking advantage of this, in the rest 
of this section we consider only a scalar position, with the understanding 
that the results can be easily expanded to the other coordinates for higher¬ 
dimensional cases. 

We assume the agents are communicating with some of their peers. In 
the most general case, the communication channels would be unidirectional, 
and therefore a directed graph is used to model the communication structure. 
Each vertex represents an agent and each edge a communication channel. 
The adjacency matrix of this graph, Ar = [aik] £ 1^”; is a binary matrix. 
The aik element is 1 if agent i receives information from agent k and 0 if 
not. The diagonal elements, an are taken as 0 too. The f-th row sum of Ar 
is called the in-degree of the vertex i [20], and it is represented by dt. The 
in-degree is also the number of peers from which the i-th member receives 
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information. They are organized in a set called the informers of agent i, 
which is denoted by J\fi. With the number of informers of each agent, the 
in degree matrix of the communication topology is created. It consists of 
a diagonal matrix that contains the in-degree of node i in its i-th diagonal 
term: A = diag (5i, (52,..., 5n) 

It is assumed that the communication process is affected by two time 
delays. The first delay is in the position communication channels and the 
second one in the velocity information exchange. The different time delays 
in the velocity and position channels is justified considering that different 
sensors are used to measure these variables, and therefore different sensing 
delays will be added on top of the delays inherent to the communication 
channels. Since we are considering that the members of the group are all 
alike, assuming the delays are uniform across the network and constant also 
makes sense. It is also assumed that each agent is able to process its own 
information fast enough to consider the so called self delay negligible. This 
means, each agent has its current information about position and velocity, 
but receives only delayed information from its peers. This assumptions are 
based on the characteristics of the physical platform currently available to 
perform a practical implementation of the methodology described in this 
paper. 

The control law is defined in this form: each agent takes the delayed 
position information it receives from its informers, computes the average and 
compares its own current position to that average; this error is penalized 
with a proportional control gain P G M^. A similar task is performed with 
the velocity information received, using a derivative control gain D G R'*'. 
Mathematically, the control law is represented by: 

•‘■(. 6,=^-'4 

for f = 1, 2 ..., n. 

Using a concatenation of the positions and velocities of the agents as the 
state vector z = [xiXiX 2 X 2 ■■■ XnXnY'i the dynamics of the multi-agent 


( 1 ) 
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system, under the control law ([T]), can be presented in state space as: 


z = 



+ (^C(8) 

+ fc® 


0 

-P 

0 0 
P 0 

0 0 
0 D 




n) 



T2) 


( 2 ) 


where 1^ is the identity matrix of size n x n, 0 represents the Kronecker 
product operation and the matrix C € is the product of the inverse of 

the in-degree matrix and the adjacency matrix: C = A“^Ar. 

The complexity of the dynamics in ([2]) increases rapidly with the number 
of agents in the group, quickly becoming intractable. In order to simplify 
the stability analysis with respect to the time delays, the methodology in¬ 
troduced in [21], extended to the case of directed topologies, is used here. 
It is based in the following lemma. 

Lemma 1. Factorization Property. The characteristic equation of the 
system m can always be expressed as the product of a £ second and m fourth 
order factors as: 


£+m 


Q {s,P,D,Ti,T2) = n {s,P,D,Ti,T 2 ,Xi) = 
i=l 

(His^ + Ds + P-Xi {Dse-^^ ^ ")] 


\i=l 
^ l-\-m 


[s^ + 2Ds^ + {D^ + 2P) + 2DPs + P^ 


( 3 ) 


\i=£+l 


- 23? (Ai) (s^ + Ds + P) [Dse -^^" + ") + 




= 0 


where Xi represent the eigenvalues of the C matrix, which has £ real eigen¬ 
values, denoted by i = 1, 2,, £ and m complex conjugate eigenvalue pairs, 
{Xi, A*), i = £ + 1, £ + 2,... £ P m. Then n = £ + 2m. 

Proof The proof is omitted here due to space constraints. It can be 
found in [TUirni] 
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As it is shown in [22], if the communication topology is connected and 
has at least one spanning tree, the matrix C has 1 as an eigenvalue with 
multiplicity one. This property introduces a factor in ([3|) of the form: 

qi (s, P, D, Ti,T 2 ) = s^ + Ds + P- {Dse-^^ ^ *) (4) 

which always has a root at s = 0 regardless of the time delays, i.e., it 
always presents a rigid body mode. This factor is known to dictate the 
dynamics of a weighted average of the states of the agents, which we call 
the weighted centroid and dictates the value upon which the agents would 
agree if consensus is reached. Without loss of generality, this eigenvalue is 
denoted as Ai. The other i + m — 1 factors in ([3|) represent the disagreement 
dynamics: if they are all stable, the consensus is reached. These results are 
well known in the literature, and are stated here as lemmas for completeness. 
The proofs of these lemmas can be found in m- 

Lemma 2. Group behavior. Assume the communication topology has at 
least one spanning tree. Then, the agents in the group reach a consensus if 
and only if the factor m is marginally stable and all the remaining factors 
of (13|) are stable. 

Lemma 3. Topologies without spanning trees. If the given communi¬ 
cation topology does not have a spanning tree, the control logic described by 
(HD cannot result in consensus. 

3 Stability Analysis Using CTCR 

The results of the previous section show that agents reach consensus if all 
factors in m are stable. Section 0] shows that this is also a necessary con¬ 
dition to ensure the formation stability. The general form of the factors 
is: 

91 (s) + 92 {s)e~'^^ ^ -h 93 (s) g4 (s) * 

which is a multiple time delay system with commensuracy (terms like 2ri) 
and delay cross-talk (terms like ti T 2 ). In order to find the delay combi¬ 
nations that render the system stable, we deploy the Cluster Treatment of 
Characteristic Roots paradigm [l8l[23], which we describe next. 

It is well known that the characteristic roots of systems like are 
continuous functions of the delays, then the stability of this class of systems 


6 


can only change when there are roots on the imaginary axis. The delay 
combinations that create imaginary roots can be classified in two classes, 
defined as follows: 

Definition 1. Kernel curves pq: The curves that consist of all the points 
exhaustively, which cause an imaginary root s = joj and 
satisfy the constraint 0 < r^o; < 2tt are called the kernel hypercurves. The 
points on these curves contain the smallest possible delay values that create 
the given imaginary root at the frequency oj. 

Definition 2. Offspring Curves p: The curves obtained from the kernel 
curves by the following pointwise nonlinear transformation: 




are called the offspring curves. 

From a point on the kernel or offspring curves, the toots can move to the 
left or right half of the complex plane, making the system stable or unstable. 
This direction of crossing is defined as the root tendency: 

Definition 3. Root Tendency The root tendency indicates the direction 
of transition of the imaginary root (to the right or to the left half of the 
complex plane) as only one of the delays increases by e, Q < e < 1, while the 
other delay remains constant: 




Root tendencies are 1 for stabilizing and +1 for destabilizing root crossings 
across the imaginary axis. 

The propositions upon which CTCR is based are the following: 

Proposition 1. Small number of kernel hypercurves: The number of 
kernel hypercurves is manageable small: for a LTI-TDS of state dimension 
n, that number is hounded by 

Proposition 2. Invariant root tendency property: Take an imaginary 
characteristic root, ju, caused by any one of the infinitely many grid points 
on the kernel and offspring hypercurves in {ti,T 2 ) € domain defined by 
expression ([6]) . The root tendency of these imaginary roots remains invari¬ 
ant so long as the grid points on different offspring hypercurves are obtained 
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by keeping one of the delays fixed. That is, the root tendency with respect 
to the variations of Ti is invariant from the kernel to the corresponding off¬ 
spring as the other delay Tk is fixed. 

The proofs of these propositions can be found in [1811^ 1^ . 

To deploy CTCR, the first step is to find the kernel curves with their 
corresponding root tendencies. This is followed by the generation of the 
offspring curves. All these curves partition the domain of the time delays 
in regions of stability and instability. To find whether a certain region 
corresponds to stable behavior, we start finding how many unstable roots 
the non delayed system has. The delays are then increased, and every time a 
kernel or offspring curve is crossed, the number of unstable roots is increased 
(if RT = 1 for that curve) or decreased (if RT = — 1) by two. Those regions 
having zero unstable roots, of course, are declared as the stable operating 
regions. 

The question now is how to find the kernel curves. For this, several 
different methods can be used, like the Kronecker sum [23], the Rekasius 
substitution [T8|, the Spectral Delay Space (SDS) based on Rekasius substi¬ 
tution [19] or the SDS based on half-angle tangent substitution [9], which 
was used in this work. Due to space limitations, we do not include details 
of the approach we followed here, and we refer the reader to PfTO] where it 
is explained. 

4 Formation Control 

The consensus protocol in ([T]) is designed to bring the state of the agents 
to a common value. If a position state is considered, that would imply a 
collision. We can transform this behavior into a formation control structure 
by just adding a constant forcing term to the consensus protocol in ([1]). In 
order to design this term, we must consider the state transformation that 
decouples the system and factorizes its characteristic equation, as stated in 
Lemma [T] and detailed in m- The matrix used in this transformation, T, 
is the one that takes C to its Jordan canonical form. Under the transfor¬ 
mation 4 = (T~^ (g) I 2 ) z, the dynamics ([2]) is decoupled into i second order 
subsystems and m fourth order subsystems. The characteristic equations 
of these subsystems are the factors of the global characteristic equation as 
declared in (l3ll . 

The state space representation of each subsystem, including the extra 
forcing term we are introducing, has the general form of: 

(i) = Mi {t) + {t-Ti) + Bi2$i (t - T 2 ) + (Pi (8) 
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where the matrices A, Bij and B 2 i have the following form for subsystems 
corresponding to real eigenvalues of C: 


A = 


0 1 
-P -D 


Bo; — \i 


0 0 
P 0 


B2i — Aj 


0 0 
0 D 


(9) 


and the following form for those corresponding to complex conjugate eigen¬ 
value pairs: 


A= I 


Bo = ( J 


B2i — 1 Ji 


0 

-P 

0 0 
P 0 

0 0 
0 D 



( 10 ) 


with Ji being: 


IHe (Aj) - Jm (A*) 
Jm(Ai) 91e(Ai) 


( 11 ) 


The forcing terms in ([8]) have the form ipi = [0 (j)i]^ for second order factors 
and ipi = [0 (pii 0 <)) 2 i]^ for fourth order factors. These forcing terms are 
applied only to the disagreement dynamics, in order to create the spacing 
among the agents. The centroid dynamics is left undisturbed, and it dictates 
where the center of the formation will be located at the end of the exercise. 

In order to study the effect of the forcing terms in ([8|), we take the 
Laplace transform and apply the final value theorem, arriving to 


1 
0 

for second order factors and to 


r 91e(A,)- 

1 -9m (Aj) 


0 

0 

1 (ku 

9m (Ai) 

Tie (Ai) - 1 

[ <t>2i 

0 

0 


P(|A 

.|2 + i_2iHe (Ai)) 


Cioo — 


B(Ai-l) 


(12a) 


(12b) 


Notice that if we make the forcing terms equal to zero, the steady state 
disagreement is zero, leading to consensus consensus. When the forcing 
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terms are present, both equations (HID imply constant relative positions 
among the agents and zero relative velocities. 

With the values of it is possible to obtain the corresponding steady 
state values for the agents positions, using z^o = (T (8)I2)^cxd- Of course, 
the formation design process does not start from but from z^o, the de¬ 
sired locations of the agents in the formation. Then, is found from = 
(T~^ (gi I 2 ) Zoo and (fT^ is used to find the value of the forcing terms in the 
transformed domain, with which the forcing vector y? = [O 0 ^ 

is created. The actual forcing terms, applied in the non-transformed domain 
are found using F = (T (8) I 2 ) 

With the formation control term added, the dynamics ([2D hnally be¬ 
comes: 


z = 



+ 

+ (c(g) 


0 

-p 

0 0 
P 0 

0 0 
0 D 




z{t) 



n) 


^ z (t - T2) F 


(13) 


5 Feedback Linearization 

The consensus protocol of section [2] and the formation control algorithm 
of section |T| were devised for second order holonomic agents, in which the 
motion along different coordinates is independent. This would be the case 
of some types of omnidirectional mobile robots. In this work, however, 
we propose the use of non-holonomic wheeled robots, more specihcally, of 
differential drive mobile robots, DDMR, like the iRobot Create [25]. In 
order to apply the protocol of previous sections to this class of robots, some 
transformation must be undertaken hrst. 

The kinematics of a DDMR is similar to that of the unicycle. The 
variables used for its description are the coordinates of the center and the 
orientation of the wheel, pictured in Fig. [1] Since the wheel is restricted 
to rolling without slipping, the relationship between the speed, v, and the 
angular rate w of the unicycle and the derivatives of its state variables is 
given by: 


X 


cos 6 

0 ■ 



1- 

1_ 


sin 9 

0 

0 

1 


V 

Ul 
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Figure 1: Variables used to describe the unicycle dynamics. 


It is well known that (jl4|) does not satisfy Brockett’s condition [26], and 
therefore there is no continuous, linear time invariant feedback law that can 
stabilize the system. The consensus protocol ([I|) cannot be directly applied. 
However, using a feedback linearization technique, proposed by m , we can 
find an equivalent linear system in which the second derivative of the state 
equals the control input, as was assumed in Section [2j 

To start the feedback linearization of m, consider a new state vector 
w = [u)i, W 2 ]^■ This vector is selected such that the derivatives of its 
components are equal to the derivatives of x and y. 

Wl 

W2 

Taking the second derivative of w we get: 

ivi 
W2 

If we make the second derivative of w equal to the system’s input, we obtain 
the equivalent linear system: 


V cos 9 — VLO sin 9 

V sin 9 VLO cos 9 


(16) 



X 



y 



V cos 9 
vsin9 


= V 


cos 9 
sin 9 


(15) 


ivi 
W2 

By using (fTHl) and (fT7|) . we obtain a relation between the input to the linear 
system and that of the non-linear system: 


Ul 

U2 


(17) 


Ul 


cos 9 

—usin0 


V 

U2 


sin 9 

V cos 9 


UJ 
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Inverting the transformation (IlSp we obtain the input that must be applied 
to the robot for it to behave as a linear system: 


V 


UJ 



cos 9 sin 9 
- i sin 9 i cos 9 



Ul 


^2 _ 


(19) 


Using (fT8]l and (USD, system da is equivalent to m , i.e., X = wi, X = wi, 
y = W2 y y = W2. 

Care must be taken when using (|19l) . since it is not defined for v = 0. 
This is known to be an structural problem of dynamic extensions of wheeled 
mobile robots m- It should be taken into account when designing control 
laws for the equivalent linear system. 


6 Simulation Examples 

In this section we consider 6 agents interacting under the communication 
topology of Fig. [21 The C matrix corresponding to this topology has the 
eigenvalues {1, 0.46, —0.5, —0.64, —0.16±j0.21} which create four second 
order factors and one fourth order factor, according to ([3|). The number 
of agents was selected because this is the number of actual experimental 
platforms that are available for a follow-on experimental study. The com¬ 
munication topology was selected such as its C matrix presented both real 
and complex eigenvalues without repetitions. For a larger number of agents 
the analysis methodology remains the same, but the computational load 
may increase, as larger communication topologies may create more diverse 
sets of eigenvalue^]. 

By deploying the CTCR paradigm into the factors, considering control 
gains P = 1 and D = 0.5, we obtain the stability chart of Fig. [3| in the 
domain of the time delays. Each one of the factors creates some kernel, 
depicted in red, and offspring curves, depicted in blue. At each point on 
those curves the system has a pair of imaginary characteristic roots. By 
increasing either of the delays at one of those points the number of unstable 
characteristic roots changes. The shaded region in that plot corresponds to 
combinations of delays ( ti , T2) for which the number of unstable roots is 
zero. That is, this is the stable operating region. For any delay combination 
within the gray area, no factor of the characteristic equation of the system 
has an unstable root. 


full communication topology, however, is very easy to analyze, as it always creates 
only two factors, regardless of the number of agents. This fact comes from results from 
graph theory [20] and is shown in detail in |28| . 
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Figure 2: Communication topology for the example cases. 



Figure 3: Stability chart for six agents interacting under the communication 
topology of Fig. [2] and protocol ([I]) when P = 1 and D = 0.5. The shaded 
region represents stable delay combinations. 
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Figure 4: Desired formation structure. 


We wish to highlight that this chart is exact and exhaustive: no delay 
combination outside of the gray area results in stable operation, and those 
in the boundary of the gray area create marginally stable operation, with 
the agents oscillating around a fixed position. 

To test the formation control algorithm, we propose the formation struc¬ 
ture of Fig. m For this configuration, the desired .^oo values are [0.71, 1.83, 
1.41, —2.08, —0.81, 5.74] in the horizontal direction and [—1.12, —0.44, 1.41, 
1.50, —1.40, 3.72] in the vertical direction. In the simulation, the agents 
start from random positions and orientations, and all have very small (10~^) 
initial speeds to avoid the mentioned division by zero in (|19h at u = 0. Using 
the delay combination (ti, t-}) = (0.3, 0.2), which corresponds to point a in 
Fig. m the agents settle into the formation after about 30 seconds. Fig¬ 
ure [5] shows the positions in x and y direction as a function of time, whereas 
Fig. 0 shows how the agents reach the desired formation as final conhgura- 
tion. The small blue lines represent the instantaneous heading of the agents 
at the moment of the snapshot. 

An increment of the delays generally worsens the performance of the 
system. When both delays are increased to (ti, T2) = (2, 1), corresponding 
to point b in Fig. [3l the agents take more than 100 seconds to reach the 
final state, as can be observed in Fig. [Tj 

Increasing the delays, however, is not always bad. An interesting, and 
rather counter-intuitive, fact that Fig. [3]shows is that, in some regions of the 
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Figure 5: Positions of the agents, in x and y direction as a function of time 
for a delay combination corresponding to point a in Fig. [3l 
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Figure 6: Traces of the motion of the agents for a delay combination corre¬ 
sponding to point a in Fig. [3l 
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Figure 7: Positions of the agents, in x and y direction as a function of time 
for a delay combination corresponding to point b in Fig. [3l 

delay space, stability can be recovered by increasing one, or both, delays. 
For example, the delay selection (ri, T 2 ) = (1, 2), marked by point c in 
Fig. El is in the unstable region. By either increasing T 2 to 5.5 seconds 
(point d) or ri to 3.5 seconds (point e), the stability is recovered. This is 
a very important result, since it provides an extra degree of flexibility to 
the control designer. Reducing the delays is normally not possible without 
changing the hardware, but increasing them is simple. 

This idea of delay scheduling was hrst discussed in [SUED]. It basically 
consists on a skillful increase of the delays to maintain stability or improve 
the performance of a system. In the case under study, we would be in¬ 
terested in the delay combination that produces the fastest convergence to 
the formation. For this, we need to calculate the dominant root of the sys¬ 
tem, and its associated time constant, for every delay combination within 
the stable operating zone of Fig. El This can be done only with numerical 
procedures, like the QPmR algorithm [31] . The stability analysis using the 
CTCR paradigm is a necessary and crucial first step for delay scheduling, 
as it is the only methodology that exhaustively declares the stable regions 
of the system, and therefore provides a finite area in which to search for the 
smallest time constant. 

For the system under study. Fig. El shows how the real part of the dom¬ 
inant root of the system changes as function of the delays within a domain 
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Figure 8: Real part of the dominant root of the system as a function of the 
delays. The black plane corresponds to 5ft (s) = 0. 


that includes most part of the stability region. The boundary depicted in 
Fig. [3] can also be seen in Fig. [8] as the intersection of the surface with the 
plane defined by 5ft (s) = 0. 

The deepest valleys of the surface in Fig. [8] correspond to the smallest 
time constants, and therefore to the fastest formation convergence rates. 
The plot shows that the best delay combinations are those in the vicinity of 
the origin, and in a small region for ri G (0, 1) and ri G (4, 6). This fact is 
better observed in Fig. [9l which represents a contour map of the surface in 

Fig. El 

With the information contained in Fig. El and Fig. [9l the designer can 
properly tune the delays to achieve a desired response time for the system. 
This is the core idea of the delay scheduling. 

7 Conclusions 

This paper considered the problem of cooperative formation control for non- 
holonomic agents in the presence of time delays. We deployed a formation 
control algorithm based on a linear consensus protocol. We consider two 
time delays: one in the exchange of position information and the other in 
the exchange of velocity information. Since the control law was originally 
devised for linear, holonomic systems, a feedback linearization technique is 
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Figure 9: Contour plot of the surface in Fig. [HI 

used to adapt the control logic to the class of agents under study. 

The stability of the formation control depends on the stability of the 
consensus protocol. In order to find the delay combinations that lead to 
stable operation, the Cluster Treatment of Characteristic Roots (CTCR) 
paradigm was used, before using CTCR, a state decomposition that leads 
to a factorization of the characteristic equation of the system is performed 
to simplify the problem. With this combination of factorization plus CTCR, 
it is possible to obtain, in a very efficient manner, an exact and exhaustive 
stability map with respect to the delays. A simulation shows that the con¬ 
vergence to formation is guaranteed in the regions declared by the proposed 
methodology. 

We also discussed briefly the concept of delay scheduling and how the 
delays can be used to tune the time the agents take to reach the desired 
formation. 

The codes used to create the examples can be requested via email to the 
authors or can be downloaded from https://db.tt/JhOOzQCf. 

Future work includes the practical implementation of these algorithms 
in a team of mobile robots. For this, the Process Control and Automation 
Laboratory of the Universidad Santo Tomas is working on the development 
of an experimental platform based on the iRobot Create development kit [25] . 
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